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TECHNICAL FIELD 

[0001] The present invention relates to a method and system for generating 
random numbers with a specified probability density function that is faster and more 
efficient than current methods and system. 

BACKGROUND OF THE INVENTION 

[0002] Methods for generating random numbers with a specified probability 
density function (PDF) are well known. One such method is the Monte Carlo Method. 
The Monte Carlo method provides approximate solutions to a variety of mathematical 
problems by performing statistical sampling experiments on a computer. The method 
applies to problems with no probabilistic content as well as to those with inherent 
probabilistic structure. 

[0003] While the conventional Monte Carlo method is often an acceptable 
sampling method, it is typically accurate only after it has converged in statistics. This 
convergence can take a very large number of random samples to provide the desired 
accuracy. For example, a simulator utilizing the Monte Carolo method may include 
many loops (for example, "For Loops" would be used in C), nested or otherwise, to 
generate an accurate statistical sample with multiple random variables. Running such 
loops with the many iterations can take hours or days, depending on the speed of the 
computer or processor. Further, the amount of time may increase more than linearly 
with the number of samples or the number of loops required to achieve convergence. 
Effective random-number generation is therefore desirable. 
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[0004] Other prior techniques for generating random numbers with a specified 
PDF utilized "brute force. 11 For example, a random number sequence uniformly 
distributed between 0 and 1 was first generated by calling on known system-provided 
routines. One difficulty with this approach is that it can take a very large amount of 
data before a sufficiently uniform distribution is realized depending upon application 
requirements. In cases where non-uniformly distributed random numbers are required, 
a controlled mapping was typically implemented in a conventional way to "shape" the 
numbers with the desired PDF. This could be done, for example, by building a lookup 
table, again with a brute-force mapping. This again can be extremely time consuming. 

[0005] Thus, it would be desirable to provide a random number generation 
method that is quicker than prior methods and also reduces the number of samples that 
must be generated to achieve the desired statistical accuracy. 

SUMMARY OF THE INVENTION 

[0006] It is therefore an object of the present invention to provide a random 
number generation method and system that requires less time than prior methods. 

[0007] It is a further object of the present invention to provide a random 
number method that reduces the number of samples needed to obtain a given overall 
accuracy for low-occurrence events. 

[0008] It is still another object of the present invention to provide a random 
number generation system that is both systematic and efficient. 

[0009] In accordance with the above and other objects of the present invention, 
a random number generation having a cumulative density function (CDF) is provided. 
(The CDF is the integral of the PDF). According to the method, a set of discrete 



numbers between an upper limit and a lower limit is generated. The set of discrete 
numbers are uniformly spaced between the upper limit and the lower limit. Each 
individual number from the set of random numbers is mapped to a corresponding 
value on a cumulative density function curve to provide a number of samples. The 
samples are then scrambled to randomize the order of the samples. 

[0010] As an additional step, the number of samples can be subjected to 
companding in certain applications. The step of companding helps provide a 
sufficient representation of a low-probability area, while decreasing the total number 
of samples. 

[0011] These and other features and advantages of the present invention will 
become apparent from the following description of the invention, when viewed in 
accordance with the accompanying drawings and appended claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0012] Figure 1 is a graph schematically illustrating the mapping of uniform 
random numbers to generate a given cumulative density function in accordance with a 
preferred embodiment of the present invention; 

[0013] Figure 2 is a graph schematically illustrating the identification of a low- 
probability area of the CDF for companding in accordance with a preferred 
embodiment of the present invention; 

[0014] Figure 3 is a flow chart schematically illustrating the sequence and 
operation o f a r andom n umber g eneration s ystem a nd m ethod i n a ccordance w ith a 
preferred embodiment of the present invention; 
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[0015] Figure 4(a) is a graph illustrating a sample without companding; 

[0016] Figure 4(b) is a graph illustrating a sample with companding; and 

[0017] Figure 5 is a flow chart schematically illustrating the sequence and 
operation of an exemplary application of a random number generation system and 
method in accordance with a preferred embodiment of the present invention. 

BRIEF DESCRIPTION OF THE PREFERRED EMBODIMENT 

[0018] The present invention is intended to allow random numbers of any 
distribution to be generated, such as through Monte Carlo simulation, in a significantly 
reduced period of time as compared to prior systems and methods. The disclosed 
system may be used in a lab environment with general purpose computers or may 
alternatively be implemented in commercial products with special-purpose processors. 
Moreover, the present invention may be utilized for a variety of applications for which 
random number generation is or can be used. 

[0019] Referring now to the Figures and specifically Figure 3, the preferred 
system and method 10 is preferably comprised of three components. The first 
component is an initial sample generation component 1 2, which generates numbers 
with immediate uniform distribution. In accordance with the number generation 
component 12, a set of numbers {R} is generated. The set of random numbers is 

preferably generated such that they have a uniform distribution between 0 and 1. 

~'f~- 

Equivalently, instead of a uniform distribution between 0 and 1, random percentages 
may be generated to have a uniform distribution between 0% and 100%. 

[0020] The first component of the invention, as generally indicated by 
reference number 12, improves distribution of random numbers by generating evenly 
spaced values for {Y} between 0 and 1. This process provides a uniform 
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representation across the interval of [0, 1]. It gives a smooth appearance to 
performance curves compiled from processing random numbers. With conventional 
generation methods, it would take a much larger number of samples to achieve the 
same curve smoothness. Alternatively, instead of generating random uniform samples 
for {R}, R can be replaced with an evenly spaced set of values {Y}. This provides a 
curve with an even smoother appearance. 

[0021] In the second component of the invention, shown by reference number 
14 of Figure 3, thereafter, a set of values {X} are generated in ascending order with a 
specified PDF and a corresponding set of descending CDF values. This is 
accomplished through a lookup table listing X vs. Y values on the CDF curve, an 
illustrative exampie o f which is shown in Figure i. The designator X is a random 
number with a CDF Y, and Y(X) is the probability that a randomly selected number 
from the set will be greater than or equal to X. In order to look-up the X value for a 
given input value R, in {R}, the conventional concept is to compare the randomly 
quantified number Rj with members of {Y} to find the closest discrete value Y k to Rj 
and therefore the corresponding value Xr. In other words, a value in the random set 
source set {R} is selected and mapped or looked-up on the CDF curve/table to 
determine the corresponding value. This process is repeated with each member of {R} 
to get a set of numbers with the specified CDF. 

[0022] As shown in Figure 1, the CDF starts with a numerical value of 1 .0 and 
decreases to 0. (For our convenience and without loss of generality, this ascending 
order of the CDF is reversed from the normal convention, which has a CDF beginning 
at 0 and increasing to a numerical value of 1). By arranging {R} in a monotonic 
order, as automatically provided by component 12, the process of comparing Rj to 
each sequential member of {Y}, beginning at Yj, can stop at the first value Yk which 



is found to be less than Rj. This process would then be repeated for every member in 
{R}. 

[0023] The preferred method provides a more efficient way to look up the X 
value for each input element in {R}. Initially, all values of {R} are arranged in, say, a 
descending order. Since both {Y} and {R} are in descending order, it would not be 
necessary to begin at Y\ each time to find the nearest value in {Y}. Since the next 
value of R, Rk+i, is smaller than Rr, Yk+i is smaller or equal to Yk. In fact, Yk+i is 
very close to Y K , and the search becomes very fast. Since significantly few 
comparisons with candidate numbers are required, significant overall savings in 
processing c an be achieved through utilization of the disclosed method and system. 
After the lookup has been completed, the order of {X} vaiues are scrambled, as 
generally indicated by reference number 16. Order randomization may be achieved in 
one o f s everal ways. For example, it may be done by generating a pseudo-random 
number (PRN) sequence {P}, which is then truncated to the length of {X} if 
necessary. {X} is then reordered by the elements of {P} as the latter are generated 
one at a time. 

[0024] The first and second components of the invention may be combined as 
follows. First, {R} is generated by uniform stepping between 1 and 0 in descending 
order. The {X} values are then generated in ascending order with table lookup. The 
lookup process proceeds without the need for pre-sorting {R} . As set forth above, at 
the end of lookup, the {X} values are then scrambled. Again, scrambling is 
systematic and is more efficient than the brute force comparison, discussed above. 

[0025] The third component of the invention, as generally indicated by 
reference number 18, is the companding component. In many applications, the 
companding improves the representation of certain regions of interest on the PDF. 
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This concept is particularly useful for PDFs where a section of interest has a much 
lower probability. This means that very few points are generated in this region by the 
conventional method, which therefore requires a much larger number of experiment 
samples to guarantee a sufficient representation of the low probability section. As the 
required members of samples increased, the run-time may increase even more rapidly. 
At the same time, more than enough samples are generated for the high probability 
areas of the PDF and are thus "wasted." Thus, the companding component increases 
the representation of low-probability samples without increasing the representation 
high-probability samples. 

[0026] For example, a given PDF might have 90% of its events happening 
between A and B in Figure 2, but oniy iOvo between B and C. Suppose ihc impact of 
interest is from the scarce region between B and C. For each sample generated 
between B and C as desired, 9 samples between A and B would normally be generated 
as well. This is where the method of companding comes in. With companding, 10 
times as many random points {R} w ould b e g enerated b etween 0 and 0.1 ofY,as 
would normally be generated. This is accompanied with a corresponding reduction in 
the number of samples generated between the popular region between 0.1 and 1.0 of Y 
for a given member of total samples. After all the points are generated and 
statistically processed, the over-represented section is re-adjusted by dividing the 
probability of that section by the companding factor N. Therefore, the low-probability 
samples are compressed while the high-probability samples are expanded. In effect, by 
reducing the spacing for low probability samples, more samples are present in the low- 
probability area. By increasing the spacing for the high probability area, fewer 
samples are present in the high-probability area. 
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[0027] Figures 4(a) and 4(b) illustrate two exemplary graphs (with CDF slopes 
in the conventional sense). Figure 4(a) illustrates a sample generation without 
companding (Y sampling uniformly spaced). Figure 4(b) illustrates a sample with 
companding. As shown by comparison of the Figures, the distribution of the samples 
provides a more accurate representation of the low-probability curve. 

[0028] The combination of these three components provide a very efficient 
and accurate method and system for generating a set of random number with a given 
PDF. 

[0029] Turn n ow t o F igure 5 , w hich i llustrates an exemplary application for 
the disclosed method and system. Preferred applications for the present invention 
include an impact analysis of rain and interference on satellite service availability. 
More specifically, this application allows a more thorough statistical link budget 
analysis when multiple statistical parameters must be analyzed. This has not 
previously been performed due to the huge amount of simulation that would be 
required to simulate interactions of the various parameters. By making such 
simulations feasible, the disclosed method can allow a more accurate link budget and 
provides better knowledge about a given system's capabilities, all with manageable 
computer processing time. By obtaining the total system performance statistics, worst 
case assumptions can be minimized. This will allow a more balanced design for 
optimal system performance within the limits of all resources. 

[0030] Referring specifically to Figure 5, rain and interference have been 
typically constant in past analysis, but can now be time varying. Thus, to determine 
the total impact of rain and interference, they can be combined in a Monte Carlo 
fashion. Initially, the parameter values from the link budget are obtained, as generally 
indicated by reference number 30. Thereafter, lookup tables for rain and interference 
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impact CDFs are built, as generally indicated by reference number 32, based on a 
variety of factors, as will be appreciated by one of skill in the art. 

[0031] Companded rain and/or an interference samples can then be generated 
based on the look-up tables, as generally indicated by reference number 34. 
Degradation histograms (PDFs) from rain with and without interference fire then 
formed, as generally indicated by reference number 36. The histograms for these two 
cases are normalized for companding, as generally indicated by reference number 38. 
The threshold of clear sky (CS) margin is then determined at 40. Finally, the change 
in availability due to interference can then be determined, as generally indicated by 
reference number 42. 

[0032] Having now fully described the invention, it will be apparent to one of 
ordinary skill in the art that many changes and modifications can be made thereto 
without departing from the spirit or scope of the invention as set forth herein. 



